Exact microscopic analysis of a thermal Brownian motor 
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We study a genuine Brownian motor by hard disk molecular dynamics and calculate analytically 
its properties, including its drift speed and thermal conductivity, from microscopic theory. 
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It is (believed to be) impossible to systematically rec- 
tify thermal fluctuations in a system at equilibrium. Such 
a perpetuum mobile of the second kind, also referred to 
as a Maxwell demon 0, would violate the second law of 
thermodynamics and would, from the point of view of 
statistical mechanics, be in contradiction with the prop- 
erty of detailed balance. Yet, it may require quite subtle 
arguments to explain in detail on specific models why 
rectification fails. Apart from the academic and peda- 
gogical interest of the question, the study of small scale 
systems is motivated by rapidly increasing capabilities in 
nanotechnology and by the huge interest in small scale 
biological systems. Furthermore, when operating under 
nonequilibrium conditions, as is the case in living organ- 
isms, the rectification of thermal fluctuations becomes 
possible. This mechanism, also referred to as a Brown- 
ian motor 0], could furnish the engine that drives and 
controls the activity on a small scale. In this letter we 
propose a breakthrough in the theoretical and numerical 
study of a small scale thermal engine. Our starting point 
is the observation that one of the basic and most popu- 
lar models, namely the Smoluchowski-Feynman ratchet 
0,0,01, is needlessly complicated, and can be replaced 
by a simplified construction involving exclusively hard 
core interactions. Its properties, including speed, diffu- 
sion coefficient and heat conductivity, can be measured 
very accurately by hard disk molecular dynamics and can 
be calculated exactly from microscopic theory. 

In Fig. we have schematically depicted the con- 
struction originally introduced by Smoluchowski in his 
discussion of Maxwell demons and re-introduced with two 
compartments at different temperatures by Feynman . 
One compartment contains a ratchet with a pawl and 
a spring, mimicking the rectifier device that is used in 
clockworks of all kinds. The macroscopic mode of op- 
eration of such an object generates the impression that 
only clockwise rotations can take place, suggesting that 
this construction can be used as a rectifier of the im- 
pulses generated by the impacts of the particles in the 
other compartment on the blades with which the ratchet 
is rigidly linked. As Feynman has argued, such a recti- 
fication is only possible when the temperature in both 
compartments is different. We now introduce a model 
which is at the same time a simplification and general- 



ization of this construction. First, we can dispose of the 
pawl and spring in the ratchet and consider any rigid 
but asymmetric object. An example with a cone-shaped 
object in one compartment and a flat "blade" or "sail" 
located in the other one is illustrated in Fig. Sec- 
ond, we replace the single rotational degree of freedom 
with a single translational degree of freedom (Figs. ^ 
and^). Third, we restrict ourselves to two-dimensional 
systems. Finally, the substrate particles in the various 
compartments are modeled by hard disks which undergo 
perfectly elastic collisions with each other while their cen- 
ter collides elastically with the edges of the motor [^. 

We first report on the results obtained from molecular 
dynamics for two different realizations of our motor. The 
first one, referred to as arrow/bar or AB is inspired by 
the above discussion. It consists of one triangular-shaped 
arrow in the first compartment and a flat bar in the other 
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FIG. 1: (a) Schematic representation of the Smoluchowski- 
Feynman ratchet, (b) Similar construction without a pawl 
and spring, (c) Two dimensional analogue referred to as the 
AB motor. The motor is constrained to move along the hori- 
zontal x-direction (without rotation or vertical displacement). 
The host gas consists of hard disks whose centers collide elas- 
tically with the engine parts. The shape of the arrow is de- 
termined by the apex angle 29o and the vertical cross section 
S. Periodic boundary conditions are used in the computer 
simulations, (d) A symmetric construction referred to as the 
A A motor. 
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FIG. 2; Probability density P(x, t) for the position x of mo- 
tor AB at times t = 1000 (shaded) and t = 4000 (open). 
Inset left: mean square displacement versus time. Inset right: 
velocity correlation function. 



(Fig. nj;). The other motor, called arrow/arrow or AA, 
consists of an identical triangular-shaped arrow in both 
compartments (Fig. QJi). Both units of the motor are 
constrained to move as a rigid whole along the horizon- 
tal x-direction as a result of their collision with the hard 
disks in the two compartments. The initial state of the 
hard disk gases corresponds to uniform (number) densi- 
ties pi and p2 and Maxwellian speeds at temperatures 
Ti and T2, in the compartments 1 and 2 respectively, 
{Ub = 1 by choice of units). The boundary conditions are 
periodic both left and right, and top and bottom. Un- 
less mentioned otherwise, the following parameter values 
are used: Each compartment is a 1200 by 300 rectangle 
and contains 800 hard disks (mass to = 1, diameter 1), 
i.e., particle densities pi = P2 = 0.00222. Initial temper- 
atures were set to Ti = 1.9 and T2 = 0.1. The motor 
has a mass M = 20, apex angle 29o — 7r/18, and vertical 
cross section S = 1. The averages are taken over 1000 
runs. 

When the temperatures are the same in both com- 
partments, Ti = T2, no rectification takes place. In 
fact. Fig. 13 shows that the motor undergoes plain Brow- 
nian motion, with average zero speed, exponentially de- 
caying velocity correlations and linearly increasing mean 
square displacement. The corresponding friction and dif- 
fusion coefficient obey the Einstein relation. On the other 
hand, as soon as the temperatures are no longer equal, 
the motor spontaneously develops an average systematic 
drift along the x-axis. The amplitude and direction of 
the speed depend in an intricate way on the parame- 
ters of the problem. In particular, the average speed in- 
creases with the temperature difference and the degree of 
asymmetry (decreasing ^o) and decreases with increasing 
mass of the motor roughly as 1/M (see the lower panel 
in Fig. OJ. Note furthermore that the observed average 
speed can be very large, i.e. comparable to the thermal 
speed yJkBT/M of the motor. The AA motor has a pe- 
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FIG. 3: Upper panel: Position of the motor as a function of 
time. The thin solid curve shows a typical trajectory. All 
other curves represent the average {x{t)). The thick solid line 
is the equilibrium case (Ti — T2 = 1) for motor AB. The 
dotted and dashed curves correspond to the nonequilibrium 
situation for respectively motor AA and AB (Ti = 1.9, T2 — 
0.1). The situation with switched temperatures for AB is the 
dashed curve with a negative velocity. Lower panel: Average 
velocity of motor AB as a function of its mass M. Inset: 
average speed (avg. over 2000 runs) of motors AA and AB as 
a function of the initial temperature difference Ti — Tb (Ti -I- 
T2 = 2 fixed) . The theoretical results O and predict lower 
speeds than the simulations. However, when the magnitude 
is scaled, the theoretical curves (dotted lines) fit well with the 
simulations. 



culiar behavior, resulting from the fact that both units 
are identical. Whereas equilibrium is usually a point of 
flux reversal, the velocity now displays a parabolic curve 
as a function of Ti — T2 with a minimum equal to zero 
at the equilibrium state Ti = T2. It is clear from its 
symmetric construction that, at least when pi = P2, an 
interchange of Ti with T2 can not modify the speed so 
that the latter has to be an even function of Ti — T2, cf. 
Fig. El 

We finally note that the observed systematic speed 
does not persist forever. Indeed, the motion of the motor 
along the x direction is a single degree of freedom that 
allows for (microscopic) energy transfers hence thermal 
contact between the compartments, a fact that was over- 
looked by Feynman in his analysis and first pointed out 
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FIG. 4: Exponential decay of the temperatures to a final com- 
mon value (Tfinai = (Ti + T2)/2 = 1.) in motor AB and con- 
comitant disappearance of the average drift speed. To enlarge 
the conductivity, a small mass M = 1, a large vertical cross 
section S — 10, and a large apex angle 6o = 7r/9 are used. 



in |3, Ig]- As a result one observes that the tempera- 
tures in both compartments converge exponentially to a 
common final temperature with a concomitant reduction 
and eventual disappearance of the systematic motion as 
shown in Fig. 0] While this feature has already been 
documented in detail in other constructions 0], we have 
focused here on conditions in which this conductivity is 
small and the compartments sufficiently large so that one 
reaches a quasi-steady state with a well defined and mea- 
surable average drift velocity. 



To obtain analytic results from microscopic theory, 
which are asymptotically exact, we focus on the situation 
in which the compartments are infinitely large while the 
densities of the hard disk gases are extremely low (more 
precisely, the so-called high Knudsen number regime re- 
quires that the mean free path is much larger than the 
linear dimensions of the motor units). In this limit, each 
compartment, characterized by its particle density pi and 
temperature Ti, acts as an ideal thermal reservoir. We 
will furthermore assume that all the constituting units 
of the motor are closed and convex. Under these cir- 
cumstances, the motor never undergoes recoUisions and 
the assumption of molecular chaos becomes exact [ir)| . 
The probability density P{V, t) for the speed V = (F, 0) 
of the motor thus obeys the following Boltzmann-Master 
equation: 

dtP{V,t) - 

j dr [W{V - r, r)P{V - r,t) ^ W{V, r)P{V, t)] (1) 

W{V,r) is the transition probability per unit time for 
the motor to change speed from V to V + r due to the 
collisions with the gas particles in various compartments. 
The explicit expression for W{V, r) follows from elemen- 
tary arguments familiar from kinetic theory of gases, tak- 
ing into account the statistics of the perfectly elastic col- 
lisions of the motor, constrained to move along the x- 
direction, with the impinging particles: 



de 



dVr 



dvyp,(t>^{v^,Vy)UF,{e){V-v)-e{B)H[{V-vye{9)]5 r + —B{9)iV - + VyCotO) 

) IVl 

m 



Here, the sum over i runs over all the different 
compartments, Li is the total circumference of the 
ith unit of the motor, B{e) = 2Msm^e/{M + 
ms\v? 9)^ H[x\ is the Heaviside function, and e[9) — 
(sin 0, — cos 9) is the unit vector normal to a surface 
at angle 9, 9 £ [0,27r], the angles being measured 
counter-clockwise from the horizontal axis. (t>i{vx,Vy) — 
me-xp[—m{vl + Vy)/2kBTi]/2nkBTi is the Maxwellian 
velocity distribution in compartment i. The shape of 
any closed convex unit of the motor is defined by the 
(normalized) probability density F{9) such that F{9)d9 
is the fraction of its outer surface that has an orientation 
between 9 and 9 + d9. Note that (sin 6*) — {cos9) — 0, 
where the average is with respect to F{9), a property 
resulting from the requirement that the object is closed. 

The Boltzmann-Master equations ^ and Q can now 
be solved by a perturbation expansion in y^m/M, fol- 



lowing a procedure similar to the one used in one- 
dimensional problems s uch as the Rayleigh piston 
or the adiabatic piston 12]. The details are somewhat 
involved and will be given elsewhere _13j. To lowest or- 
der in the perturbation, the Master equation reduces 
to a Fokker-Planck equation equivalent to the following 
linear Langevin equation: 



(3) 



with rji independent Gaussian white noises of unit 
strength and 
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the friction coefficient experienced by the motor due to 
its presence in compartment i. We conclude that at 
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this order of the perturbation the contributions from 
the separate compartments add up and are each - taken 
separately - of the linear equilibrium form. In par- 
ticular the motor has no (steady state) drift velocity, 
{V) — 0. It does however conduct heat. In the case 
of two compartments 1 and 2, the heat flow per unit 
time between them is, as anticipated in Ref. |3, given 
by a Fourier law: Qi^2 = ^(ri — T2) with conductivity 
K = kBjij2/[M{ji +72)]- One also concludes from Q 
and that the (steady state) velocity distribution of the 
motor is Maxwellian, but at the effective temperature 
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At the next order of perturbation in -y/ m/M, the cor- 
responding Langevin equation becomes nonlinear in V 
while at the same time the Gaussian nature of the white 
noise is lost, a feature well-known from the Van Kam- 
pen expansion The most relevant observation 
is the appearance, at the steady state, of a non-zero drift 
velocity: 
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In agreement with previous arguments, the A A motor, cf. 
(|7| . always moves in the same direction, namely the di- 
rection of the arrow. Furthermore it is an example where 
one can increase the asymmetry to generate a maximum 
drift speed. The limit | (sin'^ 9) \ — (sin^ 9) is reached here 
when 9o ^ 0, which corresponds to an infinitely elon- 
gated and sharp arrow in both compartments. Due to 
strong finite size effects (e.g. sound waves among oth- 
ers), the agreement between the theoretical results O 
and ^ and the computer simulations is only qualitative: 
the theory predicts speeds which are typically 20 — 40% 
lower. However, Eqs.Q and © can be fitted to the sim- 
ulation results by appropriately rescaling the magnitude 
of the velocity (see Fig. OJ, indicating that their depen- 
dencies on the parameters, M, Ti and T2, are in good 
agreement with the simulations. 

In conclusion, we have provided a detailed analytic 
and numerical study of a simplified version of the 
Smoluchowski-Feynman ratchet, including an exorcism 
- based on microscopic theory - of its operation as a 
Maxwell demon. 

This work was supported in part by the National Sci- 
ence Foundation under Grant No. DMS-0079478. 



This speed is of the order of the thermal speed of the 
motor, times the expansion parameter -y/m/Af , and fur- 
ther multiplied by a factor that depends on the de- 
tails of the construction. Note that the Brownian mo- 
tor ceases to function in the absence of a temperature 
difference (T^ = Tcff,Vz) and in the macroscopic limit 
M ^ 00 {{V) ~ 1/-M). Note also that the speed is scale- 
independent, i.e., independent of the actual size of the 
motor units: {V) is invariant under the rescaling Li to 
CLi. To isolate more clearly the effect of the asymmetry 
of the motor on its speed, we focus on the case where 
the units have the same shape in all compartments, i.e. 
Fi{9) = F{9). In this case T^s is independent of F{9) 
and the drift velocity is proportional to (sin'^ 9) / {siv? 9), 
with the average defined with respect to F{9). The latter 
ratio is in absolute value always smaller than 1, a value 
that can be reached for "strongly" asymmetric objects as 
will be shown below on a particular example. 

We now turn to a comparison between theory and sim- 
ulations. From the general result ©, one obtains the fol- 
lowing expressions for the speed of the two motors that 
were studied: 
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